ON THE CENTERS OF THE WEIGHT HOMOGENEOUS 
POLYNOMIAL VECTOR FIELDS ON THE PLANE 



Jaume Llibre 

Universitat Autonoma de Barcelona, Departament de Matematiques 
08193 Bellaterra, Barcelona, Catalonia, Spain 
jllibre@mat.uab. cat 

Claudio Pessoa 

Universidade Federal de Uberlandia, Faculdade de Matematica 
Campus Santa Monica - Bloco IF, Sala 1F118 
Av. Joao Naves de Avila, 2121, 38.408-100, Uberlandia, MG, Brasil 
pessoaOfamat . ufu.br 

Abstract. We classify all the centers of a planar weight-homogeneous poly- 
nomial vector field of weight degree 1, 2, 3 and 4. 

1. Introduction and statement of the main results 

In the qualitative theory of real planar polynomial differential systems two of the 
main problems are the determination of limit cycles and the center-focus problem; 
i.e. to distinguish when a singular point is either a focus or a center. The notion 
of center goes back at least to Poincare in [T7j. He defined it for a vector field on 
the real plane. This paper deals with the classification of the centers for a class of 
polynomial differential systems. The classification of the centers of the polynomial 
differential systems started with the quadratic ones with the works of Dulac [7j, 
Kapteyn [10], [11], Bautin [4], Zoladek [19], ... see Schlomiuk [18] for an update on 
the quadratic centers. There are many partial results for the centers of polynomial 
differential systems of degree larger than 2, but we are very far to obtain a complete 
classification of all centers for the polynomial differential systems of degree > 3. 
We deal with polynomial differential systems of the form 

x = P(x,y), 

y = Q{x,y), 

where P and Q are polynomials in the variables x and y with real coefficients. 

We say that system (JTJ) is weight-homogeneous if there exist s — (sx,S2) £ N 2 
and d € N such that for arbitrary AeM+ = {AeK:A>0}, 

P(\ s ix,\ s *y) = \ s i- 1+d P(x,y), Q(X^x,X S2 y) = X S2 - 1+d Q(x,y). 

We call s — (si, S2) the weight exponent of system (TTJ) and d the weight degree with 
respect to the weight exponent s. In the particular case that s = (1,1) systems ([T]) 
are exactly the homogeneous polynomial differential systems of degree d. 
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A singular point p of system (fT]) is a center if there is a neighborhood of p fulfilled 
of periodic orbits with the unique exception of p. The period annulus of a center is 
the region fulfilled by all the periodic orbits surrounding the center. We say that a 
center located at the origin is global if its period annulus is M 2 \ {(0, 0)}. 

Let U be an open subset of M 2 . Here a nonconstant analytic function H : U — > K 
is called a first integral of system ([1]) on U if it is constant on all solutions curves 
(x(t),y(t)) of the vector field X associated to system ([!) on U; i.e. H(x(t), y(t)) = 
constant for all values of t for which the solution (x(t),y(t)) is defined in U. Clearly 
H is a first integral of the vector field X on U if and only if 

dH dH 
XH = P—+Q— = 
ox ay 

on U. 

The main goal of this paper is to classify all centers of the weight-homogeneous 
planar polynomial differential systems (fTJ) of weight degree 1, 2, 3 and 4 with P 
and Q coprime. 

For doing that we characterize first the normal forms of all the weight-homogeneous 
planar polynomial differential systems of weight degree 1, 2, 3 and 4. Thus from the 
definition of weight-homogeneous polynomial differential systems ([1} with weight 
degree 1 , 2, 3 and 4, the exponents u and v of any monomial x u y v of P and Q are 
such that they satisfy respectively the relations 

SlU + S2V = Si, SiU + S 2 V = S 2 , 

for weight degree f ; 

silt + s 2 v = si + 1, siu + s 2 v = s 2 + 1, 

for weight degree 2; 

silt + s 2 v = si + 2, siu + s 2 v = s 2 + 2, 

for weight degree 3; and 

silt + s 2 v = si + 3, siu + s 2 v — s 2 + 3, 

for weight degree 4. Moreover taking into account that we only consider the cases 
with P and Q coprime, it is easy to check that the systems: 

(i) with weight degree f are the following ones with their corresponding values 
of s: 



(2) s = (l,l) 



x = a 10 x + a 01 y, 
y = b w x + boiy, 



(3) s=(l,p): X . ~ a h WX ' PA _ h withpeN, p>l, 

y = b p0 x p + b iy, 

(ii) with weight degree 2 are the following ones with their corresponding values 
of s: 



(4) s = (l,l) 

(5) a = (1,2) 



x = a 20 x 2 + a n xy + a 02 y 2 , 



y = b 20 x 2 + bnxy + b 02 y 2 

x = a 20 x 2 + a iy, 
y = b 30 x 3 + buxy, 



(G) 
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(iii) with weight degree 3 are the systems with weight degree 2 and additionally 
the following ones with their corresponding values of s (see [5]): 



(7) 8= (1,1) 

(8) s = (l,2) 

(9) s = (l,3) 



x = a 30 x 3 + a 21 x 2 y + ai 2 xy 2 + a 03 y 3 , 

y = b 30 x 3 + b 2 ix 2 y + b 12 xy 2 + b 03 y 3 , 

x = a 3 QX 3 + CLuxy, 

y = b 4Q x 4 + b 21 x 2 y + b 02 y 2 , 

x = a 30 x 3 + a 01 y, 

y = b 50 x 5 + b 21 x 2 y, 



(iv) with weight degree 4 are the systems (j4|), ([5]) with weight degree 2 and 
additionally the following ones with their corresponding values of s: 



(10) 


s = 


(1,1): 


X 

V 


= C140X 
= b^x* 


+ a 31 x 3 y - 
+ b 31 x 3 y-\ 


f a 22 x 2 y 2 - 
- b 22 x 2 y 2 ^ 


f ai 3 xy 3 + a 4y 
- h 3 xy 3 + b 04 y 4 , 


(11) 


s = 


(1,2): 


X 

y 


= O4o X 4 

= b 50 x 5 


+ a 21 x 2 y - 
+ b 3l x 3 y^ 


\- ao 2 y 2 , 
- b 12 xy 2 , 




(12) 


s = 


(1,3): 


X 

y 


= (140X 4 

= beox 6 


+ auxy. 

+ b 31 x 3 y + b 02 y 2 , 




(13) 


s = 


(1,4): 


X 

y 


= CI40X 

= b 70 x 7 


+ amy, 
+ h\x 3 y, 






(14) 


s = 


(2,3): 


X 

y 


= auxy, 

= b 30 x 3 + b Q2 y 2 , 






(15) 


s = 


(2,5): 


X 

i) 


= aoiy, 
= b 40 x 4 . 









For studing the centers of systems (i), (ii), (iii) and (iv), we first simplify the 
weight-homogeneous planar polynomial differential systems of weight degree 3 with 
weight exponent (1,1) having 8 parameters to some normal forms with at most 4 
independent parameters. After using these normal forms we characterize which of 
these systems have a center. The proposition that gives us this normal forms is the 
following (see [6] for the proof). 

Proposition 1. For any cubic homogeneous system there exists some linear trans- 
formation and a rescaling of the independent variable which transforms the system 
into one and only one of the following canonical forms: 



(16) * = Pix 3 +p 2 x 2 y+p 3 xy 2 +f,y 3 ^ 
y = (ix° + pix^y + p 2 xy z + p 3 y°, 

. s x = pix 3 + [p 2 + 3a)x 2 y + (p 3 - 6a)xy 2 - 6ay 3 , 

(17) . 2 , / o \ 2 , / , n \ 3 with a = ±1; 
y = Pix y + {p2- 3a)xy + (p 3 + 6a)y°, 

x = pix 3 +p 2 x 2 y + (p 3 + 2)xy 2 - 4y 3 , 

y = pix 2 y +p 2 xy 2 + (p 3 - 2)y 3 , 
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(22) 



(19) 



(20) 



(21) 




( 23 n x = xipxx 2 +p 2 xy +p 3 y 2 ), 

y = y(pix 2 y + P2xy+p 3 y 2 )- 

In what follows instead of working with system ([7]) we shall work with the equiv- 
alent normal form ([16)) - (f2"3"]) that have atmost four parameters. 
The main results of this paper is the following one. 

Theorem 2. The following statements hold. 

(a) Systems ^ has a center at the origin if and only ifboi — — aio an d Qio^oi - 
aoifrio > 0- System ([3]) has no centers. 

(b) Systems Q and ([6]) have no centers. System ([5]) has a center at the origin 
if and only if a ib 30 ^ and (2a 2 Q + &n) 2 + 8(6 30 a i - &na 2 o) < 0. 

(c) Systems ©, CES), (HB, OS), (HU), d20j and ([23]) have no centers. System 
l^fy has a center at the origin if and only if 001650 7^ 0, (3a3o + 621) 2 + 
12(650001 — 621030) < and 3a3o + 621 = 0. Systems (|21[) and (|2"2"]) have a 
center if and only if p 3 = —pi- 

(d) Systems (PH]), ((III), 02]), Gil and © have no centers. System §T3§ has a 
center at the origin if and only 1/001670 7^ and (4a4o + 631) 2 + 16(670001 — 
631^40) < 0. 

The centers of the planar weight-homogeneous polynomial vector fields of weight 
degree 1 are characterized by statement (a), of weight degree 2 by statement (6), 
of weight degree 3 by statement (c) and of weight degree 4 by statement (d) . 



In this section we assume that P and Q are coprime homogeneous polynomials 
of degree m. We want to characterize the homogeneous systems (JXJ) which have a 
center at the origin. Since P and Q are coprime it follows that the origin is the 
unique (real and finite) singular point of system ([1]) . 

System (p} in polar coordinates (p, 9) defined by x = pcos9, y = psin^, becomes 



that the function g(6) is either positive, or negative for all 6 £ S 1 ; otherwise if 
g(6*) = for some 9* E S , then the straight line 9 = 9* is invariant by the flow 
of system ([T]), in contradiction with the fact that the origin is a center. So the 
homogeneous polynomial xQ(x,y) — yP(x,y) of degree m + 1 has no real factors 
of degree 1, and consequently m must be odd. In the rest of this section we only 
consider homogeneous systems ([1]) having a center at the origin. So the function 
g(9) is either positive, or negative for all 9 & S 1 . Hence system ((T|) is equivalent to 



2. Homogeneous centers 



(24) 
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the equation 



dp _ f(0) 



d9 g{9) 



P- 



Separating the variables p and 6 of this equation, and integrating between and 9 
we get that any solution of the previous differential equation is 



From (|25|) it follows immediately the next proposition. 

Proposition 3. Suppose that P and Q are coprime homogeneous polynomials of 
degree m. Then the origin of system ^ is a global center if and only if the polyno- 
mial xQ(x, y) — yP(x, y) has no a real factors of degree 1 (in particular m is odd) 
and 



In this section we assume that P and Q are coprime polynomials and let m = 
maxjdegP, degQ}. Then we say that system ([T]) has degree m. 

Consider the polynomial vector field X = (P, Q) associated to system (TTJ) . The 
study of the local phase portrait at the singular points of the vector field X is a 
problem almost completely solved. In most of the cases one can know which is 
the behavior of the solutions in a neighborhood of a singular point. The only case 
that remains open is the monodromic one. In this case the orbits turn around the 
singular point. The difficulty is in distinguishing when the orbits spiral toward or 
backward the singular point (i.e. when the origin is a focus) and when the origin 
is a center. This problem is know as the center-focus problem, see |15j for a survey 
on this problem. 

Suppose that the origin is a singular point of X. Let j(t) be an orbit of X 
defined in a neighborhood of the origin that tends to it when t tends to +oo and 
such that linit^+oo 7(i)/||7(t)|| € S 1 , where S 1 is the unit circle, with || • || the 
Euclidean norm. In this case 7 is said a characteristic orbit. We may also consider 
orbits tending to the origin as t tends to —00 but we can always change the sign 
of the parameter t in system (fTJ and assume that t tends to +00. We have that 
the origin is a monodromic singular point of X if there is no characteristic orbit 
associated to it. 



We write P(x,y) = Y17=n P i( x ' v)> Q( x >v) = 2£L n Qi( x > v)> where Pi and Qi 



are the respective homogeneous part of degree i of P and Q, with 1 < n < m. It 
is possible that P n (x,y) or Q n (x,y) is null, but both of them cannot be null. A 
characteristic direction for the origin of X is a root u> E S 1 of the homogeneous 
polynomial xQ n (x,y) — yP n (x,y), which can be written as u> = (cos 8, sin 6) with 



It is obvious that, unless xQ n (x, y) — yP n {x, y) = 0, the number of characteristic 
directions for the origin of X is less than or equal to n+1. The following well-known 
result relates characteristic orbits with characteristic directions of singular points. 



(25) 





This proposition is well known, see for instance [6j or [12] . 



3. Monodromic singularities 



0e [0,2tt). 
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Proposition 4. Let "/(t) be a characteristic orbit for the origin of X associated to 
system ([1} and lo — lim^+oo ipj^j] ■ Then to is a characteristic direction for X . 

This proposition is proved in [3j. The reciprocal is not true, see [15] or [8]. 
In polar coordinates x = pcosO, y = psmO, and doing a change of time t i— ► r 
such that dt/dr = p n_1 , system (|T|) becomes 

f2fi s P = p(i?„(0)+pi? n+1 (0) + ---+p m -"i? m (p)), 

1 j = F n (6)+pF n+1 (6) + --- + p m - n F m (p)), 

where 

i?j (0) = Pj (cos 0, sin 9) cos + (cos 9, sin 0) sin 9, 
Fj (9) — Qj (cos 9, sin 0) cos 9 — Pj (cos 0, sin 9) sin 0, 

are homogeneous trigonometric polynomials of degree j + 1, for n < j < m. 

Let 7(i) = (p(i), #(£)) be an orbit for system fll} written in polar coordinates. We 
have that 7(t) tends to the origin if p(t) is not identically zero and lim t ^ +co p(t) = 0. 
Moreover j(t) tends to the origin spirally if it tends to the origin and linit^+oo 9{t) = 
±oo. Therefore j(t) is a characteristic orbit for the origin of X if it tends to the 
origin and lim t ^ +00 9(t) = 9* < oo. In this case we say that 9* is the tangent at 
the origin of j(t). 

The polynomial xQ n (x, y) — yP n (x, y), considered over S , in polar coordinates is 
F n (9) given in ([26]) . Hence if ^(t) = (p(t), 9{t)) is a characteristic orbit for the origin 
of X written in polar coordinates with tangent 9* at the origin, then F n (9*) = 0. 

Now we include here a well-known property of monodromic singular points, see 
[8] for an easy proof. 

Proposition 5. // the origin is monodromic, then F n (9) ^ and 9(p,9) is of 
definite sign for any (p, 9) verifying < p < e, for a certain e > sufficiently 
small. 

When the eigenvalues of the matrix DX(p) at the singular point p are complex 
and not real, we know that the origin is monodromic. Such class of singular points 
are called of linear type. If their real part is different from zero then the singular 
point is a focus, while if their real part is zero the singular point may be a center 
or a focus. This last case is the classical Lyapunov-Poincare center problem. The 
study of this problem for a concrete family of differential equations passes through 
the calculation of the so-called Lyapunov constants. The problem of how many 
constants are needed to distinguish between a center of a focus is in general open. 

When the matrix DX{p) has its two eigenvalues equal to zero but the matrix is 
not identically null, it is said that p is a nilpotent singular point. The monodromy 
problem in this case was solved in [2] and the center problem has been studied in 
[16j and [9]. Nevertheless, in general, given a polynomial system with a nilpotent 
monodromic singular point it is not an easy task to know if it is a focus or a center. 

Suppose that X has a singular point p such that the matrix DX(p) is a nilpotent 
matrix. If the singularity p is a center then p is called a nilpotent center. In this 
case using suitable coordinates system (p} can be written as 

( 27 n x = y + P 2 (x,y), 

y = Q2{x,y), 

where P2, Q2 are polynomials of degree at least 2. 

We state a theorem, proved in [2] (see also [3] and [1]), which solves the mon- 
odromic problem for nilpotent singular points. 
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Theorem 6. Consider system (|2T[) and assume that the origin is an isolated sin- 
gularity. Define the functions 

f[x) = Q 2 (x,F(x)) = ax a + 0(x a+1 ), 

<p{x) = dW(y + P 2 (x,y),Q 2 {x,y)) \ y=F{x} = bx? + 0{xP +1 ), 

where a ^ Q, a > 2, b ^ and (3 > 1, or 4>{x) = 0. Here the function y — F{x) 
is the solution of y + P2(x,y) = passing through (0,0). Then the origin of (|27p 
is monodromic if and only if a is negative, a is an odd number (a = 2n — 1), and 
one of the following three conditions holds: [3>n—l;f3 = n—l and b 2 + 4an < 0; 
and = 0. 

Once we know how to distinguish between monodromic and not monodromic 
nilpotcnt singular points, we wish to solve the center-focus problem for the mon- 
odromic ones. Our approach to this problem passes through the computation of the 
Taylor expansion of the return map near the singular point. In the next section we 
define the generalized Lyapunov constants introduced to solve the stability problem. 
Instead of working with system (|27|) we shall use the special normal form given in 
the next statement, see [3] or [1] for the proof. 

Lemma 7. An analytic vector field having an isolated nilpotent monodromic sin- 
gularity can be written as 

/ 2g ^ x = y(-l + P(x,y)), 

V = f(x) + 0(x)y + Q(x,y)y 2 , 



where all the above functions are analytic at the origin and satisfy P(0, 0) = 0, 
f(x) = X 2 ™" 1 + 0(x 2n ) for some n € N, n > 2, and <p(x) = or tj>(x) = bx* 3 + 
0(x^ +1 ) for some (3 £ N satisfying (3 > n — 1. Moreover, if (3 = n — 1 then 
b 2 -An< 0. 

4. Generalized Lyapunov constants 

Using the notation of previous section we will introduce some generalized polar 
coordinates. As a starting point recall that given any natural number n € N, the 
generalized trigonometric functions x{9) = Cs((9), y(8) = Sn(0), are defined as the 
unique solution of the Cauchy problem 

dx 

(29) 

v _ de ' 

with initial conditions x(0) = 1, y(0) = 0. 

Some properties of these functions are stated in the next proposition, see |14j . 

Proposition 8. Let n6N and let Cs{6) and Sn(0) be the generalized trigonometric 
functions determined by system ()29|1 . Then the following statement hold. 

(a) Cs 2n (6») +nSn 2 (6») = 1. 

(b) Cs(0) and Sn(9) are T -periodic functions where 

T = 2^ T{ ^ ) 



nT (2±i' 



2n ) 

(c) Sn(0 + T/2) = Sn(9), Cs(6» + T/2) = -Cs(6), Sn(-6» + T/2) = Sn(6») ano 
Cs(-9 + T/2) = -Cs(9). 



8 



JAUME LLIBRE AND CLAUDIO PESSOA 



Here as is usual T denotes the gamma function, i.e. T(z) — e~ t t z ~ 1 dt. 
We define for any given n £ N the generalized polar coordinates r and 9 of the 
real plane (x, y) £ K 2 as 

(30) x = rCs(9), y = r n Sn(d), 

where the function Sn and Cs are given through (|29|) . Notice that x 2n + ny 2n = r 2n . 
Furthermore the following equalities hold 

x 2n ~ 1 x + yy 

r = 

r 2n~l ' 

(31) 

• xy - nyx 

Introducing the coordinates ([30]) . it follows from ([31]) that system (|28| can be 
reduced to the following differential equation 

(32) de= R{r > 9)= i + H (r,ey 

where R is analytic at the origin. Thus the above equation can be written in the 
neighborhood of the origin as 

, oo 

(33) i=E^K> 

i—l 

where the functions Ri for i > 1 are T-periodic, being T the period of the general- 
ized trigonometric functions. Now we can define the generalized Lyapunov constants 
see [14]. 

Consider the solution r(0,p) of such that r(0, p) — p. It can be written as 

oo 

(34) r(9,p)=Y i u i (9)p i , 

i=l 

where Ux(0) = 1 and Wfc(0) = for all k > 2. Hence the return map is given by the 
series 

oo 

P(p)=r(T,p)=J2ui(T)p i - 

i=l 

For a fixed system the only significant term is the first that makes the return map 
differ from the identity map, and this will determine the stability of the origin. On 
the other hand if we consider a family of systems depending on parameters, each 
of the Ui(T) depends on these parameters. Thus the stability of the origin is given 
by the first V k = u k (T) ^ with u x (T) = 1, u 2 (T) = ■■■ = u fe _ x (T) = 0, which it 
is called fcth generalized Lyapunov constant. Note that system (|27p has a center if 
and only if Vi = 1 e V k = for all k > 2. 
In [T] it is proved the following result. 

Proposition 9. Consider system (|28|) under the assumptions of Lemma^ Then 
the origin is a monodromic singular point and its first generalized Lyapunov con- 
stants is 

when (3 = n — 1 and n is odd; otherwise V\ = 1 . 
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5. Weight-homogeneous centers of Weight degree 1, 2, 3 and 4 

In this section we shall prove the main results of this paper. We remark that for 
a weight-homogeneous polynomial vector field X = (P, Q) with P and Q coprime 
the origin is the unique (real and finite) singularity of X. Moreover when it is a 
center it is a global center (see |12|). 

Proof of Theorem fj| (a) . Consider system ([2]) . This system is linear and we know 
that the origin is a center if and only if 6 i = — a 10 and aio&oi — aoi^io > 0. Now 
system © does not have a center because the straight line x = is invariant by its 
flow. □ 

Proof of Theorem^ (6). Consider the vector field X = (P, Q) associated to system 
(fJJ. This vector field does not have centers by Proposition [3l because m is even. 

System (j6j does not have a center, because the origin is always a cusp, by a 
result that can be find in [3] page 362. Note that this system has the first integral 
H(x,y) = a 01 y 2 /2 - b 20 x 3 /3. 

Now we will study system ([5]). Note that 001630 7^ 0, otherwise the straight 
lines x = or y = should be invariant by the flow of the system, and so the 
origin cannot be a center. Therefore rescaling the independent variable system (|5|) 
becomes 

id/ i "20 2 
x = y + P 2 {x,y) = y^ x, 

aoi 

n 1 \ ^30 3 011 

V = Q 2 (x,y) = — x H xy. 

aoi aoi 

We have that F(x) — — a2ox 2 /aoi is the solution of y + Pi{x, y) — 0. Hence 

/(*) = eh(x,F(x)) = 630001 ; 6llfl2 v , 

aoi 

<f>(x) = div(y + P 2 (x,y),Q 2 (x,y)) \ y=F ( w) = 2a20 + bll a;- 

a i 

Thus by Theorem [S] the origin of system ([5]) is monodromic if and only if 
(2a 2 o + foil) 2 + 8(6 30 a i - 611020) < 0. 

In this case the origin is a center of ([5|), because the system above is reversible, i.e. 
it is invariant by the change of variables (a;, y, t) 1— > (—a;, y, —t). □ 

Proof of Theorem^ (c). We begin studying system l[5]). This system does not have 
centers, because the straight line x = is invariant by its flow. 

Consider system Note that 001650 ^ 0, otherwise the straight lines x = or 
y = should be invariant by the flow of this system, and so the origin cannot be a 
center. Rescaling the independent variable system © becomes 

1 a 30 3 

x = y H x , 

( 35 ) . 62i 01 2 , 6 50 5 

y = = x y H x . 

aoi aoi 

In a similar way as in the study of system ([5]) in the proof of Theorem [5] (6) , we 
have that the origin of system ([9]) is monodromic if and only if 

(3a 30 + 6 2 i) 2 + 12(6 50 a i - 6 2 ia 30 ) < 0. 
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Now by the change of variables 



/«30&21 — aoi^5o\ , fl 3° / a 30^21 — «01^5o\ 1 ^ 



y = 2 y i 

V a 01 / a oi V a 

( O30&21 — 101^50^ 
X = ~\ 2 

system ([35]) can be written as 



„3 



X, 



x = -y, 

(36) . 5 i • ( \ 3a 30 + 6 2 i 2 

y = x +sign(a i) -==== x y. 

v 030O21 — 050001 

Therefore by Proposition [5] the first generalized Lyapunov constant V\ of system 
(HI is 

^ ofr,;„-„/ n - \ 3a 3 o+&2l ^„ ^ 



V\ = exp 



ny 4 n-(sign( ao ^ i ^ oaoi ) 2 j 



Hence 3o,3o + 621 = is a necessary condition for the origin be a center of system 
(HI). But this condition is also sufficient, because in this case system (|36|) has the 
first integral H(x, y) — y 2 /2 + x 6 /6. Note that the origin is a global minimum of 
the graphic of the H . 

For finishing the proof of statement (c) we must study system ([7]) . We consider 
the normal forms of system ([7]) given by Proposition [T] 

Let X = (P, Q) be the vector field associated to the systems given by Proposition 
Q] System (fl6|) does not have a center by Proposition [3j because due to the fact 
that 

xQ(x, y) - yP{x, y) = n(x -y)(x + y)(x 2 + y 2 ), 

it has the invariant straight lines x ± y = though the origin. 
In the same away system (fT7|) does not have a center, because 

xQ{x, y) - yP{x, y) = -6ay 2 (x 2 - 2xy - y 2 ). 

Now systems (TTg|) . (|T^)) . PU]) and ([^]) do not have a center, because y = is an 
invariant straight line for these systems. 
For system (P2!?1) we have that 

xQ(x, y) - yP(x, y) = a(x 2 + y 2 ) 2 . 

On the other hand for this system we get that 

W) a ' 

where / and g are given by ([24]) . Therefore, by Proposition [3J system (|22[) has a 
center if and only if p 3 — —pi . 

Analogously for system (|2"T|) we have that 

xQ(x, y) - yP(x, y) = a(x 4 + 6fix 2 y 2 + y 4 ). 

Hence xQ{x, y) — yP(x, y) = if and only if 



= ±|y|V-3/*± W-l. 
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Thus, if fi > —1/3 the polynomial x 4 + 6px 2 y 2 + y 4 does not have a real factor of 
degree 1. On another hand for this system we get that 

^ fie) m 

g{v) V9/1 2 - la 

where / and g are given by ([24]) . Therefore, by Proposition [3] system (|2ip has a 
center if and only if p 3 = —p\ . This finish the proof of Theorem [U (c) . □ 

Proof of Theorem^ (d) . Consider the vector field X = (P, Q) associated to system 
(TTU)) . This vector field does not have centers by Proposition O because m is even. 

Systems (fT2"j) and (fTi|) does not have centers, because the straight line x = is 
invariant by its flow. 

System (fT5j) does not have a center, because the origin is always a cusp, by a 
result that can be find in [3] page 362. Note that this system has the first integral 
H(x, y) = b 4Q x 5 /5 - a 01 y 2 /2. 

Consider system (IT3"|) . Note that 001670 7^ 0, otherwise the straight lines x = 
or y — should be invariant by the flow of this system, and so the origin cannot 
be a center. Rescaling the independent variable system (|13p becomes 

. a 40 4 

x = y H x , 

( 37 ) . 631 3 & 70 7 

y = — x y H x . 

a oi a 01 

In a similar way as in the study of system ([5]) in the proof of Theorem [2] (6) , we 
have that the origin of system (|13[) is monodromic if and only if 

(4a 40 + 631) 2 + 16(6 70 aoi - 631040) < 0. 

Now if (J01&70 7^ and the origin of system (|13j) is monodromic, then it is a center 
because system ([37| is reversible, i.e. it is invariant by the change of variables 
(x,y,t) 1 ► (-x,y,-t). 

For finishing the proof we must study system (TTT]) . Note that 002650 7^ 0, other- 
wise the straight lines x = or y — should be invariant by the flow of this system, 
and so the origin cannot be a center. 

Now in polar coordinates system is given by (|26|) . where n = 2 and £2(0) = 

-a 02 sin 3 (6'). Hence, as 0(p, 8) = -o 02 sin 3 (6') + pF 3 (9) H , it follows that 6(p, 9) 

changes of sign for < p < e with e small enough. Therefore, by Proposition the 
origin cannot be a monodromic singular point of system Hence system (jlip 

has no centers. □ 
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